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The effects of short-range correlations derived from a realistic meson-exchange potential on the 
single-particle density matrix in finite nuclei are investigated by analyzing the one-body density in 
terms of the natural orbits. Basic features of these natural orbits and their spectral distributions 
are discussed. For many observables it seems to be sufficient to approximate the one-body density 
matrix in terms of those natural orbits, which exhibit the largest occupation probabilities. For 
the investigation of the high-momentum components in the single-particle density, however, it is 
important to take into account natural orbits with small occupation probabilities, originating from 
the single-particle Green function at large negative energies. 
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It is generally accepted that atomic nuclei are many-body systems in which correlations beyond the mean-field or 
Hartree-Fock picture play a significant role. Therefore it has always been a point of experimental and theoretical 
studies to find observables, which reflect these correlations in a unambiguous way. In this sense, the simultaneous 
study of density and momentum distributions in nuclei seems to be a sensitive test of nuclear models and may 
provide information on the correlations between the interacting nucleons 0]. Density and momentum distributions 
are one-body properties and therefore can be calculated from the single-particle spectral function which contains all 
the information on the single-particle structure of nuclei. Recent proton knock-out experiments employing the (e, e'p) 
reaction have shown that correlations between nucleons are responsible for a reduction of the spectroscopic strength 
. of valence hole states with respect to the uncorrelated mean-field values and for the fragmentation of the strength 
over a large excitation energy domain of the remaining nucleus §-@. 

The modern (e, e'p) performed e.g. with the MAMI facility [Q in Mainz or at NIKHEF Ej] in Amsterdam do not only 
provide information on the global spectroscopic strength but also yield results on the high-momentum components 
of the spectral function. The effects of correlations on the momentum and energy distribution of the single particle 
strength, in particular for these high momentum values, is presently a subject of discussion [p|-[Io|. Most of the 
investigations have been performed for light nuclei or infinite nuclear matter |ll|-[l7[|. Results for heavier nuclei 
are typically derived from investigations of nuclear matter assuming a local density approximation |l8-21 1. 



A question that has attracted a lot of attention concerns the possibility of describing a given momentum distribution 
by means of an independent particle model. In the IPM the nucleus is considered to be a system of nucleons moving 
without residual interaction in a mean-field or single-particle potential. In this case the nuclear wave function is 
represented by a single Slater determinant. One way to clarify this point is by studying the natural orbits (no) 
J,^2-25|. These wave functions and their associated eigenvalues (n a ) (natural occupation numbers) are obtained 



by diagonalizing the one-body density matrix for the ground state of the A particle system. Therefore the one- 
body density matrix of the correlated system can be expressed in terms of these occupation numbers n a and the 
corresponding natural orbits 0™° by 

p(r,rO = £n Q C°*MCV) (1) 

a 

As the one body density matrix is diagonal in the no representation, this basis appears to contain the most suitable 
single-particle wave functions for the calculation of the expectation value of a general one-body operator. In particular, 
if the total wave function of the system is a unique Slater determinant build with single-particle wave functions then 
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the natural orbits coincide with those wave functions and the natural occupations are equal to unity. The deviations of 
this ideal situation indicate the impossibility to obtain the one-body density matrix from a single Slater determinant. 
One can say that the effect of correlations in the calculation of one-body observables are smallest when using the no 
basis; for example, the depletion of the Fermi sea is smallest in this basis Gol. 



Performing a microscopic calculation of the one-body Green function for O we have recently demonstrated [£6 27 
that the nuclcon-nucleon correlations induced by the short-range and the tensor components of a realistic interac- 
tion |28| yield an enhancement of the momentum distribution at high momenta as compared to the Hartree-Fock 
description. This enhancement, however, is only obtained after integrating over all excitation energies of the A — 1 
particle system. On the other hand, a momentum distribution very similar to the one derived in the Hartree-Fock 
approximation is observed if one restricts the energy integration to the states of lowest energy for each set of orbital 
angular momentum I and j. Does this imply that correlation effects are visible in (e, e'p) experiments to the low-lying 
states in A ~ 1 nucleus only in terms of the spectroscopic factors, which are reduced as compared to the independent 
particle model? In order to discuss this question, we will analyze the results of the single-particle Green function 
obtained with the method described in p7| in terms of natural orbits and their occupation numbers. This analysis will 
also allow us to study the effects of correlations on the local density distribution p(r) and other one-body observables 
in 16 0. 

After this introduction, the techniques to evaluate the single-particle Green function, the density matrix and the 
natural orbits are briefly described in the next section. The discussion of our numerical results are presented in section 
3 and the main conclusions of the work are summarized in the final section 4. 

II. CALCULATION OF THE SINGLE PARTICLE DENSITY MATRIX 

The aim of our investigation is the analysis of the single-particle density matrix p mn , which is defined in a basis of 
appropriate single-particle states m, n by 

Pmn=<*o|aL«n|*0 > ( 2 ) 

where denotes the correlated ground state of the nucleus and (a n ) define the single-particle creation (annihila- 
tion) operators for the basis considered. For a nucleus like 16 O with a ground state of angular momentum J = the 
density matrix can be separated into sub-matrices of given orbital angular momentum /, total angular momentum j 
and an isospin quantum number. In our discussion we will ignore the effects of the Coulomb interaction. This means 
that the results obtained in nuclei with N = Z are identical for protons and neutrons and therefore we will drop the 
quantum number referring to the isospin. Using e.g. the momentum representation and performing a partial wave 
decomposition, the density matrix can be written as pij(k\, k 2 ), with Ij referring to the "good" quantum numbers of 
angular momentum and fci, k 2 define the absolute values of the momenta for the single-particle states m,n above. 

Within the Green function approach this density matrix can be evaluated from the imaginary part of the 
single-particle Green function by integrating 

Pij(ki,k 2 ) = ( dE-lm(gij(ki,k 2 ;E)), (3) 

J-oo 

where the energy variable E corresponds to the energy difference between the ground state of the A particle system 
and the various energies of all the states in the A — 1 system (negative E with large absolute value correspond to high 
excitation energies in the remaining system) and is the Fermi energy. 

The single-particle Green function gij or propagator is obtained from a solution of the Dyson equation 

g lj (k 1 ,k 2 ;E) =gf (k 1 ,k 2 ;E) + J dk 3 J dfc 45 j9 \k u h; E)AE tj (k 3 , k A ; E)g h (fc 4 , k 2 ; E), (4) 

where g(°' refers to a Hartree-Fock propagator and ASjj represents contributions to the real and imaginary part of 
the irreducible self-energy, which go beyond the Hartree-Fock approximation of the nucleon self-energy used to derive 
g(°\ The definition and evaluation of the Hartree-Fock contribution as well as the calculation of AE will be discussed 
below. 

The calculation of the self-energy is performed in terms of a G-matrix which is obtained as a solution of the Bethe- 
Goldstone equation for nuclear matter choosing for the bare NN interaction the One-Boson-Exchange potential B 
defined by Machleidt ( p8[ |, Tab. A. 2). The Bethe-Goldstone equation has been solved for nuclear matter with a 
Fermi momentum kp = 1.4/m . This roughly corresponds to the saturation density of nuclear matter. The starting 
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energy has been chosen to be -10 MeV. The choices for the density of nuclear matter and the starting energy are 
rather arbitrary. It turns out, however, that the calculation of the Hartree-Fock term is not very sensitive to this 
choice |30|. Furthermore, we will correct this nuclear matter approximation by calculating the 2-particle 1-hole (2plh) 
term displayed in Fig. [l]b directly for the finite system, correcting the double-counting contained in the Hartree-Fock 
term (see discussion below). 

Using vector bracket transformation coefficients pif , the G- matrix elements obtained from the Bethe-Goldstone 
equation can be transformed from the representation in coordinates of relative and center of mass momenta to the 
coordinates of single-particle momenta in the laboratory frame in which the 2-particle state would be described by 
quantum numbers such as 

\k1hj1k2hj2JT) , (5) 

where ki , ij and ji refer to momentum and angular momenta of particle i whereas J and T define the total angular 
momentum and isospin of the two-particle state. It should be noted that Eq. (|5|) represents an antisymmetrized 
2-particle state. Performing an integration over one of the ki, one obtains a 2-particle state in a mixed representation 
of one particle in a bound harmonic oscillator while the other is in a plane wave state 

/>oo 

\n1hj1k2hj2JT) = dkiklRn u h(oiki) \k1hj\k2l2j2JT) ■■ (6) 
Jo 

Here R ni ,h stands for the radial oscillator function and the oscillator length a = 1.72 fm _1 has been selected. This 
choice for the oscillator length corresponds to an oscillator energy of hui osc — 14 MeV. Therefore the oscillator functions 
are quite appropriate to describe the wave functions of the bound single-particle states in 16 O. Indeed, it turns out 
that the single-particle wave functions determined in self-consistent BHF calculations for 16 O have a large overlap 
with these oscillator functions |32). Using the nomenclature defined in Eqs. (|J) - @ our Hartree-Fock approximation 
for the self-energy is easily obtained in the momentum representation 

Sg£ 2(2 1 +l) 2 (2J + l)(2T + l)(k 1 hjin 2 l 2 j 2 JT\G\k' 1 hj 1 n 2 l 2 j 2 JT). (7) 

The summation over the oscillator quantum numbers is restricted to the states occupied in the IPM of 16 0. This 
Hartree-Fock part of the self-energy is real and does not depend on the energy. It should be noted that the approxima- 
tion, which is called Hartree-Fock in this paper is not a Hartree-Fock approach in terms of the bare NN interaction V. 
Using a realistic NN interaction V, as it is done here, such a Hartree-Fock approach would yield a very unreasonable 
description, in which the nucleons are typically unbound |33] ]. What is called Hartree-Fock approximation here, is 
defined in terms of a nuclear matter G matrix, which accounts already for NN correlations. 

In order to evaluate now the single-particle Green function and densities we consider a complete basis within a 
spherical box of a radius Rbox- This box radius should be larger than the radius of the nucleus considered. The 
calculated observables are independent of the choice of i?box, if it is chosen to be around 15 fm or larger. A complete 
and orthonormal set of regular basis functions within this box is given by 

$i«m(r) = (r\hljm) = Nuji(kir)yi jm (64>) (8) 

In this equation ^ m represent the spherical harmonics including the spin degrees of freedom and ji denote the 
spherical Bessel functions for the discrete momenta fcj which fulfill 



Using the normalization constants 



ji(hR hox ) = 0. (9) 



—, r, for/>0 



Nil = ( v^-^bo*) (1Q) 
' ^ hlvJ_ forZ = 0, 



the basis functions defined in Eq. (||) are orthogonal and normalized within the box 

rfbox 



d r (ki>l j m\r) (r\hljm) = 5u<Su>Sjj>5 m m> ■ (11) 



n 
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Note that the basis functions defined for discrete values of the momentum fej within the box differ from the plane 
wave states defined in the continuum with the corresponding momentum just by the normalization constant, which is 
■v/2/7r for the latter. This enables us to determine the matrix elements of the nucleon self-energy in the basis of Eq. 
(|J) from the results presented in Eq. ((?]). For that purpose the Hartree-Fock Hamiltonian is diagonalized 



iW 2 



E 



(ki\ ^S in + Ey \k n ) (fc„|a) y = (h\a) l3 . (12) 



Here and in the following the set of basis states in the box has been truncated by assuming an appropriate iV max . In 
the basis of Hartree-Fock states |a), the Hartree-Fock propagator is diagonal and given by 

9%\a;E)= , (13) 

where the sign in front of the infinitesimal imaginary quantity irj is positive (negative) if e^J is above (below) the 
Fermi energy. Knowing the expansion coefficients of the Hartree-Fock states in the "box-basis" (ki\a)ij it is of course 
trivial to present gfy in the basis of Eq. (||). 

The contributions to the self-energy, which are of second order in the nuclear matter G matrix are represented by the 
diagrams displayed in Figs. 0b) and [lp), referring to intermediate 2plh and 2-hole 1-particle (2hlp) states, respectively. 
These two contributions define the correction to the self-energy AS beyond the Hartree-Fock approximation in the 
Dyson Eq. (|J). These two contributions yield a complex and energy-dependent correction to the self-energy. The 
2plh contribution to the imaginary part is given by 

W^ h (h, k[; E) = E nM2 E hUj3ji Ejt I k 2 3 dk 3 J k 2 dk 4 (2J + 1)(2T + 1) 

x {k 1 lij 1 n 2 hj2JT\ G Ih^j^Uj^JT) 
x {k3hjskilijiT\ G |fc' 1 /ijin 2 / 2 j2 JT) 

E + e n2 i 2 j 2 - — 2^r) , (14) 



W5( 



where the "experimental" single-particle energies en 2 z 2 j 2 are used for the hole states (-47 MeV, -21.8 MeV, -15.7 MeV 
for S1/2, P3/2 and P1/2 states, respectively), while the energies of the particle states are given in terms of the kinetic 
energy only. The expression in Eq. ([IJ) still ignores the requirement that the intermediate particle states must be 
orthogonal to the hole states, which are occupied for the nucleus under consideration. The techniques to incorporate 
the orthogonalization of the intermediate plane wave states to the occupied hole states as discussed in detail by 
Borromeo et al. |Q have also been used here. The 2hlp contribution to the imaginary part W 2 ^ p (ki, k[; E) can be 
calculated in a similar way (see also [|4}). 

Our choice to assume pure kinetic energies for the particle states in calculating the imaginary parts of W 2plh (Eq. 
(|l4|)) and W 2hlp may not be very realistic for the excitation modes at low energy. Indeed a sizeable imaginary part 
in W 2hlp is obtained only for energies E below -40 MeV. As we are mainly interested, however, in the effects of 
short-range correlations, which lead to excitations of particle states with high momentum, the choice seems to be 
appropriate. A different approach would be required to treat the coupling to the very low-lying two-particlc-onc-holc 
and two-hole-one-particle states in an adequate way. Attempts at such a treatment can be found in Refs. pSj . 

The real parts of the 2plh and 2hlp terms in the self-energy can be calculated from the corresponding imaginary 
parts by using dispersion relations pq| . As an example we present the dispersion relation for the 2plh part, which is 
given by 



V h p h (h ,k 1 ;E) = - J ^— dE, (15) 



where P means a principal value integral. Since the Hartree-Fock contribution T, HF has been calculated in terms of 
a nuclear matter G-matrix, it already contains 2plh terms of the kind displayed in Fig. ^b). Therefore one would 
run into problems of double counting if one simply adds the real part V 2plh to the Hartree-Fock self-energy. Notice 
that T, HF does not contain any imaginary part because it is calculated with a nuclear matter G-matrix at a starting 
energy for which G is real. In order to avoid such an over counting of the particle-particle ladder terms, we subtract 
from the real part of the self-energy a correction term V c , which just contains this contribution calculated in nuclear 
matter EtII. 
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Summing up the various contributions we obtain the following expression for the correction to the Hartree-Fock 
self-energy 

AE = (V 2plh -V c + V 2hlp ) + i (W 2plh + W 2hlp ) . (16) 

Rather than calculating the various contributions to AE in the basis of plane wave states, one may also calculate 
AE in the box basis defined in Eq. (Q) or in the basis of Hartree-Fock states a, b, which we have obtained from the 
diagonalization in Eq. (|l2|). In this basis it is easy to determine the so-called reducible self-energy, originating from 
an iteration of AS, by solving 

(o| T\f {E) \b) = (o| ASy(^) \b) + <«l A ^(E) \c) g§\ci E) (c\ T\f{E) \b) (17) 

C 

and obtain the propagator, which corresponds to the solution of the Dyson Eq. ([|) from 

9lj (a, b; E) = g$\a; E) (o| T,\f{E) \b) <#>(&; E). (18) 

Using this representation of the Green function one can calculate the continuum contribution to the density matrix 
in the "box basis" from 

pfjih,^) = \ J E2hlP dElm ^<Ai|a) Ii5w (a,6;£?)<6|A a )yj , (19) 

where the integration limit E 2 hip refers to the negative threshold of 2hlp states for the Ij combination considered. 
In our studies this energy is well below the corresponding Hartree-Fock energies. Therefore the imaginary part of <?y 
in the region of integration is different from zero only due to the imaginary part in E red , which just describes the 
coupling of the quasi particle states to the 2hlp continuum. 

Besides this continuum contribution, the hole spectral function also receives contributions from the quasihole states 
p9[ . The energies and wave functions of these quasihole states can be determined by diagonalizing the Hartree-Fock 
Hamiltonian plus AE in the "box basis" 

E t^ 8m + + ASy(jB = 4 « } |fcn> (fcn|T> ^ = e ^ (fci|T> ^- • (20) 

n— 1 

Since in the present work AE only contains a sizeable imaginary part for energies E below e^ 1 , the energies of the 
quasihole states come out real and the continuum contribution is separated in energy from the quasihole contribution. 
The quasihole contribution to the density matrix is given by 

p%{k l M) = Zyi j {k x \T) l .{T\k 2 ) w (21) 

with the spectroscopic factor for the quasihole state given by |29] 

d (T| AEy(i?) |T) 



Zrij = 11- 



dE 



.qh 
'Tlj 



(22) 



Finally, the total sing le-p article density matrix is given as the sum of the continuum part of Eq. (|19|) and the quasihole 
contributions of Eq. ( J2l| ) summed over all quasihole states with an energy below the Fermi energy e p 

p l ^k 1 ,k 2 )=pt j (k 1 ,k 2 )+ E P%ikuk 2 )- (23) 

ef l} <e F 

At this stage it is now of course straightforward to determine the natural orbits and their occupation numbers by 
diagonalizing this density matrix according to Eq. ([!]) . 
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III. RESULTS AND DISCUSSION 



A. Natural orbits and occupations 



The numerical results to be discussed in this section have been obtained for each combination of orbital angular 
momentum I and total angular momentum j by using a basis of single-particle states defined according to Eqs. (||) - 
( |l0| ) for a spherical box with a radius R of 20 fm. Basis states up to A max =20 were taken into account. With these 
assumptions the density matrix for each single-particle channel Ij is represented by a matrix of dimension 7V max =20. 
The diagonalization of these matrices yields eigenvalues which are the occupation probabilities ?iy a of the natural 
orbits (see Eq.(|])) and the eigenvectors correspond to the expansion coefficients of these natural orbits {4>ija) in the 
chosen basis. 

Results for these occupation probabilities are displayed in table | ordered for each partial wave (Ij) according to 
the magnitude of the occupation probability. Here and in the following we will use the notation that the natural 
orbit identified by a=l corresponds to the orbit with largest occupation probability and also the orbits with a=2...20 
are ordered with respect to the occupation probability n a . Note that natural orbits with a larger than 4 show an 
occupation which is so close to zero, that their contribution to the total occupations for each Ij, which are also given 
in table |], is hardly visible within the accuracy used in that table. 

For those partial waves Ij, for which one orbit is completely occupied in the Hartree-Fock (HF) approach (sj/2, 
P3/2, and P1/2 in our example 16 O), one finds that the first natural orbit (a=l) exhibits an occupation probability 
very close to one, while the occupation probabilities of the natural orbits for a =2 and larger are smaller by about a 
factor of 0.01. Similar observations have been reported for the natural orbits of 3 He drops containing 70 atoms. |3^] 



It is interesting to observe that the eigenvalues for the 'occupied' orbits in |39| are in general 10-40% smaller than the 



corresponding ones calculated here for b O, indicating that correlations between ''He atoms are stronger than those 
between nucleons. For those orbits Ij which are completely unoccupied in the HF approach, like d§/ 2 , d^j 2 or /V/2 hi 
our example, the occupation probabilities are of course much smaller in agreement with the results in p9]. The ratio 
n a=2 versus n a= \ for these orbits is not as small as in the case of the orbits with occupation in the HF approach. 
As we will discussed more in detail below, this difference is due to the quasiholc contribution to the density in these 
"occupied" orbits, which is missing in the Ij channels which don't show any occupation in HF. 

The total occupation of these orbits with I > 1 is decreasing with increasing I. This could be interpreted to reflect 
the feature that orbits which have a larger excitation energy within the shell-model are less occupied. It is remarkable, 
however, that for a given I the states with j = I — 1/2 are more populated than those with j = I + 1/2, although 
the spin-orbit term of the shell-model yields less binding for the former than for the latter states. This is true for 
the example of the d states displayed in the table but also for states with I > 2, which are not given explicitly here. 
In order to understand this feature we must recall that the occupation probability for e.g. natural orbits with d 5 / 2 
symmetry in 16 O, reflects the probability to generate in knock-out reactions states of the A — 1 particle system 

\Mt 1} >^ «yl*o A) > ■ (24) 

The resulting state \^\f > is a one- hole state with respect to the correlated ground state of 16 O, it does not 
contain, however, any component with one hole in the d 5 / 2 orbit with respect to the HF ground state configuration 
for 16 O. Such configurations do not exist or do not have the particle number (A — 1)=15. The occupation probability 
of natural orbits with d 5 / 2 symmetry reflects the possibility of generating states with two- hole 1-particle or more 
complicate (n+l)-hole n-particle configurations as compared to the HF ground state configuration for 16 0. This 
probability, however, is not affected by single-particle energies of the d 5 / 2 orbits but rather determined by the number 
and energies of these (n+l)-hole n-particle configurations and the probability of these configurations to be generated 

(A) 

by applying ay to the correlated ground state 'J'q . Keeping this in mind the occupation probabilities of the natural 
orbits with I > 1 are easily understood. 

The radial shapes of natural orbits with a—1 to 3 are displayed in Fig. |[ At first sight these natural orbits look 
like solutions of a Schrodinger equation for a particle moving in a single-particle potential. They exhibit the same 
Z-dependence (4>ij a ~ r ) for small values of r, the number of nodes is essentially identical to a — 1, and their values 
decrease for large distances as it is the case for single-particle wave functions for bound states. 

In order to characterize these natural orbits with different a on a more quantitative level, we have calculated the 
square radii of these orbits 

<r 2 > lja = f r 2 drcf>f 3a (r)r 2 (25) 



G 



and show some typical results in table H. It should be emphasized that the resulting < r 2 >ij a do not increase with 
the number of nodes, i.e. the value of a, as one expects from single-particle wave functions in a potential. Assuming 
e.g. an oscillator potential the corresponding values for < r 2 > increase as a function of a proportional to (2a+l— 1/2). 
The results for the natural orbits are much less increasing with a, in the case of the I — 1 orbits the values obtained for 
a = 2 are even smaller than the corresponding results for a = 1 . This means that the natural orbits with a > 1 yield 
a density distribution much more localized in the center of the nucleus than would be obtained from corresponding 
single-particle wave functions for a local potential. 

Furthermore it is remarkable that the values of < r 2 > obtained for the a = 1 natural orbits with angular momenta 
I = 2 and I = 3 are smaller then the radii calculated for the natural orbits with large occupation probability and 1 = 1. 
This difference is due to the fact that the 1 = 1, a = 1 natural orbits are dominated by the quasihole contribution to 
the density matrix (see Eq. ^l|). The quasihole wave functions have overlaps larger than 0.99 with the corresponding 
natural orbits a = 1. The quasihole states contribute an occupation of 0.78, 0.914 and 0.898 to the total occupations 
of 0.921, 0.947 and 0.93 for the a = 1 natural orbits with lj equal to s 1 / 2 , P3/2 and P1/2, respectively. The radii of 
the quasihole states are slightly larger than the radii of the corresponding natural orbits (see table || and also figure 
^). This indicates that the continuum contribution to the natural orbits tend to have smaller radii than the quasihole 
part. Therefore it can also be understood that radii of natural orbits with I > 1, which do not contain any quasihole 
contributions, can be smaller than those for 1 = 1. 

From a more general point of view, however, one may say that for all those orbits, which are occupied within the 
naive shell-model and therefore yield a quasihole state, the radial wave function for the natural orbit with largest 
occupation probability (a = 1) is very similar to the wave function of the quasihole state. Furthermore, both wave 
functions are very close to the corresponding single-particle wave functions obtained in the HF approach (see figure 

D- . . 

Summarizing this part of the discussion we find that for a given Ij the continuum contribution to the density matrix 



(see Eq. 19) yields smaller radii than the corresponding quasihole wave function. Even if we sum the contributions 



of all natural orbits 



< r z >u= ^ a ^ (26) 

the result is smaller than just the radius of the quasihole state (see table [0]). In order to obtain the total radius of 
the nucleus one has to sum over all partial waves Ij, multiply each < r 2 >y by the total occupation of these orbits 
E Q n ija times the degeneracy (2(2j + 1)) and divide by the total particle number. Due to the depletion of states with 
I < 1 accompanied by a small occupation of orbits with I > 1 one obtains a value for the total radius (R=2.55 fm), 
which is only slightly below the result obtained in the Hartree-Fock approach (R=2.59 fm). 

These features are also displayed in figure ||, where the local density distribution obtained in HF (long dashes) 
is compared to the total density distribution (solid line). Also we present the contribution to the total density 
distribution coming from partial waves with / > 1 (short dashes). One finds that the depletion of states with small I 
(in particular / = 0) yields a reduction of the density in the center of the nucleus. The difference between the density 
distribution from the quasi-hole approximation (dashed-dottcd line) and the total result is due to the continuum 
contribution to the density. 



B. Spectral decomposition 

In the preceding discussion it has been pointed out already that the natural orbits with occupation probabilities 
close to one arc dominated by the quasihole contribution. This means that features of these natural orbits are explored 
in knock-out experiments like (e, e'p) populating the states of low energy in the daughter nucleus. In the following 
we would like to discuss the question, to which extent the analysis of experiments with large missing energies, leading 
to high excitation energies in the residual (A — 1) system, exhibits the properties of the natural orbits with smaller 
occupation probabilities. 

For that purpose we analyze the continuum contribution to the density matrix (p c , see eq.(|l^)). For the example 
of orbits with 1 = 0, the largest occupation numbers resulting from the diagonalization of this continuum part of the 
density matrix are listed in table III and the corresponding natural orbits are displayed in figure |5|. The natural orbit 



with largest occupation, obtained from this continuum part of the density matrix, yields a radial wave function, which 
is again very similar to the corresponding quasihole state. In fact, the overlap of the natural orbit displayed in the 
left part of figure ||has an overlap with the quasihole wave function of the / = state, which is as large as 0.997. This 
means that also the continuum part of the density matrix is characterized by a dominant natural orbit very similar 
to the quasihole state. 
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In a next step we split the energy integration of eq.(|T^) into 2 parts. The first interval ranges from E = —100 
MeV to the threshold of the 2h\p continuum, while the second interval covers energies from -co to E = —100 MeV. 
Restricting the energy integration to one of these two intervals we can study either the single-particle density to 
be observed in knock-out experiments with very large missing energies (interval 2) or at medium missing energies. 
Eigenvalues (i.e. occupation probabilities n a ) and eigenvectors (natural orbits) of these parts of the density matrix 
are also presented in table III and figure ^[ respectively. 



The natural orbits originating from the energy interval 2 (large missing energies) are more localized in the center of 
the nucleus than those derived from the energy interval which corresponds to lower excitation energies in the residual 
system. The overlaps of the various natural orbits with the same a, however, are still large (typically larger than 
0.95). What seems to be more significant is the relative importance of the different energy intervals for the various 
natural orbits. The occupation of the a — 1 natural orbit mainly originates from the energy interval 1, whereas the 
natural orbits with a equal to 2 or larger, collect the main part of the occupation from energies below -100 MeV in 
the integral of eq. (|l^) . 

The same features are also observed in the other partial waves as can be seen from the results displayed in table 
[iv| . In this table we give for various partial waves Ij the occupation probabilities of the natural orbits derived from 
the total density matrix, split into the contribution resulting from the quasihole part of the density matrix (denoted 
by n^ 1 ) and the continuum part (n c a ). These results show again that the a = 1 natural orbits are dominated by 
the quasihole component of the density matrix, while the other natural orbits mainly originate from the continuum 
contribution to the single-particle density. 

Furthermore we have listed for each natural orbit, 4>ij a , the mean value for the energy in the 2hlp continuum (see 
also Eq.©) 

wi ia = ^7 / ^ dEElm ( }^ (n|/ M ),. //,,■ ( k, . k->: E) ;/,- 2 |n) /; J . <27i 

yk ± ,k2 

In this equation Qij(ki,k2] E) stands for the single-particle Green function in the "box-basis" of Eqs.(|J) - (|Io| ) and 
(ki\a)ij denote the expansion coefficients of the natural orbits in this basis. 

Inspecting these mean values for the energies one finds that the absolute values of these energies increase with a, 
the index of the natural orbits chosen with respect to the occupation probability. This implies that the natural orbits, 
which have radial wave functions with a larger number of nodes (a > 1, see figure ||), "collect" the occupation at 
more negative energies than the natural orbits with a = 1. 

Also we list in table ^ the continuum contribution to the expectation values for the kinetic energies of the natural 
orbits 



J 2hlp 



= ^ / dE ( k r\ T kin\h)ij Im ( { a \kx)ij9lj{ki,k 2 ;E) (fc 2 |a) y ] . (28) 

i k 1 ,k 2 

Also the expectation values iy a increase with a. For a=l they are slightly larger than the corresponding expectation 
values for the quasihole states. This is in line with the observation made above that the continuum part of the 
density matrix tends to decrease the radius of the natural orbit with a — 1 as compared to the quasihole state. The 
expectation values for a larger than one increase very drastically with a since the corresponding natural orbits contain 
strong components with high momenta as can bee seen from the radial shapes in figure 

Concluding this part of the discussion we find that the natural orbits with a larger t i 
occupations) represent the high momentum components in the single-particle density and originate mainly from the 
imaginary part of the Green function at very negative energies. This means that the features of these natural orbits 
might be explored by analyzing knock-out experiments with large missing energies. 



lan 1 (indicating their small 



The occupation probabilities and expectation values listed in table IV can also be used to calculate the total energy 
of the nucleus 16 O via the "Koltun sum rule" 



^2*(2j + l) 



ij 



(29) 



The quasihole contribution to the total energy is calculated in terms of the quasihole energy e^y (see Eq.(|2f]|)), the 

kinetic energy for the quasihole states and the spectroscopic factors Zrij (see Eq. (p2|) ) . It is remarkable that this 
quasihole contribution yields only around 37 percent of the total energy although the quasihole part of the single- 
particle density is responsible for 89 percent of the total particle number. The largest contribution to the energy 
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in Eq. ( |29| ) originates from the continuum part of the natural orbits with a = 1 . These states contribute more than 
58 percent of the total energy although the particle number of these contributions is only around 9 percent. The 
remaining part of the total energy £, which is -3 MeV out of the total calculated E = — 81.1MeV (see also p^] ), is 
due to the natural orbits with a larger than 1. 

It is worth repeating that a large part of the density matrix can be described in terms of one natural orbit for each 
partial wave. These contributions of the a = 1 natural orbits account for 97.9 percent of the total particle number 
and for 96 percent of the total energy. Therefore the description of the single-particle density matrix in terms of these 
natural orbits with largest occupation probabilities should be sufficient for the calculation of many single-particle 
observables. It should be kept in mind that already this description in terms of one natural orbit in each partial wave 
is beyond the mean-field approximation: (i) It allows for a depletion of orbits, which are completely occupied in the 
mean-field approach, and occupation of partial waves, which are not occupied in the mean-field picture, (ii) This 
approach takes into account the most important part of the spectral distribution of the single-particle density matrix 
(see the contributions to the energy in Eq.(p9[)). The restriction to the natural orbits with a=l, however, should not 
be sufficient to obtain a reliable description of the high-momentum components in the single-particle density matrix. 
These components are described mainly in terms of the natural orbits with smaller occupation probability. 



IV. CONCLUSIONS 



Properties of the correlated wave function for the nucleus O, which includes the effects of short-range correlations 
derived from a realistic NN interaction by means of the Green function method |p7f , have been analyzed in terms 
of natural orbits. It turns out that the single-particle density matrix can very well be described in terms of a few 
natural orbits for each partial wave Ij. If one restricts the representation of the density matrix in Eq.(]l|) to only one 
orbit per partial wave (the orbit with largest occupation probability n a , here referred to as the a = 1 natural orbits), 
one obtains a density matrix, which accounts already for 97.9 percent of the total particle number. These a = 1 
natural orbits have a large overlap with the wave functions of the quasihole states as well as the Hartree-Fock (HF) 
single-particle wave functions for those orbits, which are occupied within the mean-field approximation. 

At first sight also the natural orbits in partial waves, which are unoccupied in the HF approximation, and the 
natural orbits with small occupation probabilities (a > 1) look like single-particle wave functions for a nucleons in an 
excited state of a single-particle potential. A closer inspection shows, however, that the radii for these "excited" orbits 
are much smaller than those obtained for a motion within a single-particle orbit. These a > 1 natural orbits mainly 
originate from the single-particle Green function at large negative energies (E <-100 MeV). Although the occupation of 
these orbits is much smaller than those with a=l, these orbits describe the high- momentum components of the single- 
particle density matrix. Therefore in order to study these "excited" natural orbits or the high-momentum components 
of the single-particle density, one has to analyze nucleon knock-out experiments with large missing energy. 

For the understanding of the bulk properties of nuclei, like binding energy and radius, it seems to be a good 
approximation to restrict the study to natural orbits with maximal occupation probability. It is not sufficient, 
however, to restrict the evaluation of the energy according to the "Koltun sum rule" (Eq|2^) to the contributions of 
the quasihole states. The main contributions arise from the single-particle Green function at energies deep in the 
2hlp continuum. In order to improve the present studies towards a self-consistent treatment of the single-particle 
Green function, it should be a good approximation to describe this Green function in terms of the natural orbits with 
largest occupation probabilities. 
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us (H.M.) is pleased to acknowledge the warm hospitality at the Facultat de Fisica, Universitat de Barcelona, and 
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TABLE I. Occupation of natural orbits in le O. Listed are the occupation probabilities n a (see eq.l) of those 4 natural orbits, 
which show the largest occupation probabilities in various partial waves. Also given are the total occupation probabilities, i.e. 
the sum over all n a including orbits with a > 4. 



a 


Sl/2 


P3/2 


Pl/2 


ds/2 


ds/2 


fr/2 


1 


0.9206 


0.9467 


0.9304 


0.0154 


0.0189 


0.0056 


2 


0.0127 


0.0073 


0.0083 


0.0033 


0.0045 


0.0013 


3 


0.0024 


0.0014 


0.0019 


0.0007 


0.0010 


0.0003 


4 


0.0004 


0.0003 


0.0004 


0.0001 


0.0001 


0.0000 




0.936 


0.956 


0.941 


0.0195 


00245 


0.0073 
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TABLE II. Radii of Natural Orbits. Listed are the expectation values of r 2 for natural orbits with largest occupation 
probabilities in various partial waves. Also given are the expectation value for the total single-particle density (< r 2 >y, see 
Eq. (25)) in each partial wave (normalized to 1) and the corresponding results for the Hartree-Fock (HF) and the quasihole 
states (qh). All entities listed are given in fm 2 . 



a 


«l/2 


£>3/2 


Pl/2 




d-s/2 


fr/2 


1 




4.092 


7.159 


7.535 


5.746 


5.967 


6.270 


2 




4.446 


5.761 


6.070 


7.033 


7.215 


9.070 


3 




6.978 


9.981 


10.364 


14.105 


13.563 


16.296 


4 




17.982 


20.445 


20.172 


22.776 


22.872 


26.877 


< r 2 > tj 


4.112 


7.155 


7.534 


6.370 


6.582 


7.392 


HF 




3.974 


7.435 


7.963 








qh 




4.164 


7.258 


7.653 








TABLE III. Occupation probabilities for natural orbits obtained by diagonalizing the continuum part of the density matrix 


for the I 


= partial wave. Results 


are presented restricting the energy 


integral of Eq.(19) to the intervals shown 


in the first 


line. 




















[- 


100MeV,e F ] 




[-oo, -lOOMeV] 




[-oo, e F ] 


1 






0.09390 




0.04920 




0.14118 


2 






0.00465 




0.00627 




0.01232 


3 






0.00058 




0.00121 




0.00222 


4 






0.00006 




0.00018 




0.00037 








0.09919 




0.05693 




0.15612 
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TABLE IV. Expectation of natural orbits in various partial waves lj. Listed are the part of the occupation probabilities 
arising from the quasihole (n q J l ) and the continuum (n„) contribution to the single-particle density, the mean value for the energy 
(uiija, see Eq.(27)) and the kinetic energy (tij a , Eq.(28)) in the 2hlp continuum. Furthermore we also give the single-particle 
energy (e qh ) and the kinetic energy t qh of the quasihole states, if they exist. 



a 


n qh 


< 




Jijc [MeV] 


tijc [MeV] 






S1/2 e 9 = 


-34.30 MeV, 


t q = 11.23 MeV 




1 


0.77988 


0.14068 




-87.13 


12.49 


2 


0.00006 


0.01268 




-118.22 


50.28 


3 


0.00002 


0.00234 




-126.71 


95.49 






ah 

PS/2 ^ = 


-17.90 MeV, 


t q = 18.06 MeV 




1 


0.91420 


0.03251 




-87.22 


23.70 


2 


0.00009 


0.00718 




-122.75 


64.32 


3 


0.00001 


0.00139 




-126.26 


98.74 








d'o/2 






1 




0.01539 




-95.65 


44.44 


2 




0.00332 




-115.89 


89.77 


3 




0.00073 




-128.97 


143.07 
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FIG. 1. Graphical representation of the Hartree-Fock (a), the 2-particle 1-hole (2plh, b) and the 2-hole 1-particle contribution 
(2hlp, c) to the self-energy of the nucleon 

FIG. 2. The radial shape of natural orbits with lj = S1/2, P3/2 and d 5 / 2 - Ordered with respect to the occupation number, 
results are displayed for a — 1 (solid line), 2 (long dashes) and 3 (short dashes) 

FIG. 3. Orbital shape of the natural orbit (a — 1, solid line) the Hartree-Fock (short dashes) and the quasihole wave function 
in the Si/ 2 partial wave 

FIG. 4. Density distribution of nucleons in ie O as a function of the distance r from the center of the nucleus. The result 
obtained in the HF approximation (long dashes) is compared to the one obtained within the Green function approach (solid 
line). Also displayed are the contributions of the quasi-hole states to this total density (dash-dotted line) and the density 
originating from the occupation of states with I larger than 1 (short dashes) . 

FIG. 5. Natural orbits resulting from the diagonalization of the continuum contribution to the single-particle density in I = 
partial wave. The natural orbits obtained from the energy interval [— lOOMeV, 6f] in Eq.(19) (short dashes) and those from 
the interval [—00, — lOOMeV] (long dashes) are compared to those from the whole energy range (solid line). The left part of 
the figure shows the natural orbits with a = 1, while the right part displays those for a = 3. 
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